We perform a seismic study of the young massive (3 Cephei 
star HD 203664 with the goal to constrain its interior struc- 
ture. Our study is based on a time series of 328 new Geneva 
7-colour photometric data of the star spread over 496.8 days. 
The data confirm the frequency of the dominant mode of the 
star which we refine to f\ = 6.02885 cd~ x . The mode has a 
large amplitude of 37 mmag in V and is unambiguously iden- 
tified as a dipole mode (£ — 2) from its amplitude ratios and 
non-adiabatic computations. Besides f±, we discover two ad- 
ditional new frequencies in the star with amplitudes above 4a: 
f 2 = 6.82902 cd" 1 and / 3 = 4.81543cd~ 1 or one of their 
daily aliases. The amplitudes of these two modes are only be- 
. - - , tween 3 and 4 mmag which explains why they were not de- 
■ tected before. Their amplitude ratios are too uncertain for mode 
| identification. We show that the observed oscillation spectrum 
^vq ■ of HD 203664 is compatible with standard stellar models but 
\ that we have insufficient information for asteroseismic infer- 
q ■ ences. Among the large-amplitude (3 Cephei stars, HD 203664 
stands out as the only one rotating at a significant fraction of its 
critical rotation velocity (~ 40%). 
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Abstract. 



1. Introduction 

The /3 Cephei stars are a homogeneous group of oscillating 
B0-B3 stars that have been studied as a class for more than a 
century now. Stankov & Handler (2005) recently compiled an 
overview of the observational properties of this group of stars. 
The oscillations of (3 Cephei stars are explained in terms of the 
k mechanism operating in the ionisation layer of the iron-peak 
elements (Cox et al. 1992; Kiriakidis et al. 1992; Moskalik 
& Dziembowski 1992; Dziembowksi & Pamyathnykh 1993). 
Given that mainly low-degree low-order pressure and gravity 
modes are excited, these stars are good potential targets for 
in-depth seismic studies of the interior structure of massive 
(i.e. pre-supernova) stars. Indeed, the luminosity classes of the 
known (3 Cephei stars range from V up to I (Stankov & Handler 
2005; see also Waelkens et al. 1998) and theory predicts the oc- 
currence of oscillations for this whole area in the HR diagram 
(Pamyatnykh 1999). 

Recent progress in the seismic interpretation of selected 
f3 Cephei stars was remarkable in the sense that standard stellar 
structure models are unable to explain the oscillation data for 
the best-studied stars: HD 129929 (Aerts et al. 2003), v Eridani 
(Pamyatnykh et al. 2004, Ausseloos et al. 2004), 12Lacertae 
(Handler et al. 2005, Ausseloos 2005). Pamyatnykh et al. 
(2004) have suggested to include radiative diffusion processes 
in a new generation of stellar models in an attempt to re- 
solve the excitation problem in v Eridani. This has not yet been 
achieved so far. These three well-studied (3 Cephei stars are all 
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slow rotators with vsini below 40kms 1 and a mass between 
9 and 12 M . 

In view of these recent achievements, and in an attempt to 
obtain similar results for a star with higher surface rotation ve- 
locity, we have selected one of the most rapid rotators among 
the large-amplitude (3 Cephei stars for a long-term photometric 
monitoring programme on which we report here. 

The star HD 203664 (B0.5V, m v =8.59) was discovered to 
be a new (3 Cephei star by Aerts (2000), who derived one 
frequency of 6.0289 cd" 1 from the HIPPARCOS photometry. 
This frequency was confirmed by her in 49 Geneva measure- 
ments spread over about one year taken with the P7 photometer 
attached to the 0.7m Swiss telescope at La Silla observatory. 
The scarce multicolour data set did not allow discrimination 
between i = 1 or 2 for the spherical degree of this oscillation 
mode but did seem to exclude a radial mode. HD 203664 is 
among the top ten of the class as far as photometric amplitude 
is concerned, with a value of ~ 30mmag in the Geneva V filter 
(Aerts 2000). It is by far the most rapid rotator among the large- 
amplitude members (see Fig. 4 of Stankov & Handler 2005), 
with vsini=200km s _1 derived from high-resolution spectra by 
Little et al. (1994). HD 203664 also happens to be one of the 
very few class members which is situated at high galactic lat- 
itute (I = 61.°93,6 = -27.°46) at a distance of 3200pc, 
thanks to which high-precision spectroscopic data is available 
and carefully analysed (Little et al. 1994). 

In this paper, we report the findings of our observational 
study of HD 203664 in an attempt to contribute to a better un- 
derstanding of (3 Cephei stellar structure models in a diversity 
of such type of stars. 

2. Data description and stellar parameters 

We included HD 203664 in the long-term photometric moni- 
toring programme of pulsating stars performed with the 1 .2m 
Mercator telescope at Roque de los Muchachos in La Palma, 
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Canary Islands. In this framework we obtained 328 Geneva 
7-colour high-precision photometric measurements between 
HJD 2452085.6 and HJD 2452582.3. The time span of these 
new data is 496.8 days. The integration times were typically 
4 minutes, resulting in a precision of about 7 mmag per mea- 
surement in U and 6 mmag in V. Aerts (2000) had obtained 
already 49 datapoints for the star between HJD 2450391 and 
HJD 2450790 with the same instrument but attached to the 
0.7m Swiss telescope at La Silla. The Southern and Northern 
Geneva standard star systems are carefully calibrated so that 
measurements in both hemispheres should be compatible, even 
over a long baseline. All reduced data are provided in Table l. 1 

The basic stellar parameters of HD 203664 were derived 
from different sources. Aerts (2000) used the old Geneva pho- 
tometry and positioned the star in the HR diagram with respect 
to the (3 Cephei star instability strip (see her Fig. 1). Using stan- 
dard stellar models published by Schaller et al. (1992), she thus 
derived a mass of 13. %Mq. We have recomputed the estimates 
using the same method as in Aerts (2000) from the average 
value of the 6 Geneva colours for all new data and find refined 
values of log T ofr = 4.47 ± 0.01, log g = 3.9 ± 0.3. 

It is well-known, however, that fundamental parameter es- 
timates for B stars from high-resolution spectroscopy often re- 
sult in a lower effective temperature and gravity (see, e.g. De 
Ridder et al. 2004 for a discussion about this for the (3 Cephei 
star ^Eridani). Moreover, in the case of HD 203664 we can- 
not rely on an accurate value of the parallax, as shown by the 
large discrepancy between the result of 0.32 mas by Little et al. 
(1994) and of 2.23±1.08mas from HIPPARCOS (Perryman et 
al. 1997). Little etal. (1994) derived the stellar parameters from 
high-resolution spectra and used these to estimate the distance. 
They find log T cff = 4.447, log g = 3.7 and a mass of 14 M Q 
(based on evolutionary models by Maeder & Meynet 1988), 
as well as a normal (i.e. solar) abundance for B stars in our 
vicinity. Unfortunately, these authors did not provide error es- 
timates. Finally, the few f3 Cephei stars with accurate seismic 
modelling have always ended up outside their observationally 
determined error box in effective temperature and gravity, to 
the cooler and less massive part (Thoul et al. 2003 for 16 Lac, 
Aerts et al. 2003 for HD 129929 and Pamyatnykh et al. 2004 
and Ausseloos et al. 2004 for v Eridani). 

The final estimate of the parameters we conservatively 
adopt for HD 203664, based on all these arguments, is 
log T cff = 4.45 ± 0.02 and logg = 3.8 ± 0.2, while we do 
not use any constraint at all on its luminosity. 

3. Frequency analysis 

We searched for frequencies in the Geneva U, B, V filters with 
the method outlined in Scargle (1982). The results were similar 
in the B and V filter, so we only provide the detailed analysis 
for the U and B filter. We accepted frequencies as long as their 
amplitude is more than 4 times the noise level, which corre- 
sponds to a 99.9% confidence level of having found an intrin- 
sic variation rather than a frequency due to noise (Breger et al. 

1 available only in electronic form at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr (130.79.128.5) 



1993, Kuschnig et al. 1997). This criterion is common prac- 
tise among asteroseismologists these days (e.g. Handler et al. 
2003, 2004, 2005). The noise level was computed by averag- 
ing the periodogram peaks in the range [0, 10] c d _1 after final 
prewhitening. 

The accuracy of the frequencies was calculated as 
<T f ~ a/^NAT (Home & Baliunas 1986, Montgomery & 
O'Donoghue 1999) where the proportionality constant is of or- 
der 1 depending on the author (see Cuypers 1987 for a discus- 
sion), A is the amplitude of the frequency /, N is the number 
of data points, T is the total time base and a is the average 
error on each individual measurement. We estimated the lat- 
ter by computing the standard deviation of the noise after final 
prewhitening, and found it to be 7.1 mmag for the Geneva U 
filter, 6.5 mmag for B and 7.1 mmag for V. 

During a first step, we confirmed the dominant frequency 
already found by Aerts (2000) from the HIPPARCOS photom- 
etry. This frequency appeared clearly in all three considered fil- 
ters: /i = 6.02885(2) cd -1 (where the uncertainty of the last 
digit is given in parenthesis). There was strong aliasing due to 
the single-site nature of the data (see top panels in Fig.Q but 
as the HIPPARCOS data gave us the same value without any 
daily alias we were sure that we picked the correct frequency 
for the dominant mode. 

After prewhitening, we encountered strong frequency 
peaks at ncd -1 . These peaks turned out to due to the slight 
difference in the average magnitude between the older 0.70m 
data and the 1 .2m Mercator data, which introduces daily aliases 
of the yearly periods due to the observing seasons. This small 
difference in average magnitude probably results from the fact 
that the standard stars used for the Mercator measurements are 
fainter than those used for the 0.70m Swiss telescope at La 
Silla. In order not to be disturbed by these daily aliases, we ig- 
nored the 49 older measurements in our subsequent frequency 
analyses. 

Prewhitening with fx then led to a clear second frequency 
fa = 6.82902(13) cd" 1 and its aliases in the U filter. This 
frequency and its aliases fulfilled safely the amplitude require- 
ment (Fig.Q. In the B filter, we found the yearly alias of fa: 
f 2 = 7.83176cd~ 1 . The HIPPARCOS photometry was of no 
use to help discriminate among the aliases as of this stage, be- 
cause there are no significant frequencies in that dataset after 
prewhitening with fx . We proceeded by considering fa as well 
as f 2 in a biperiodic fit together with fx, but this did not help to 
discriminate between the two options. Therefore, we concluded 
that either fa or f 2 is the true second frequency. Phase plots 
for fx and fa for the three Geneva UBV filters are provided in 
Fig.|2 The curves in the lower panel are indistinguishable in 
quality from those for f 2 . 

After subsequent prewhitening with either fa or f 2 we con- 
tinued the search for new frequencies. In this way, we found 
fa = 4.81543(14) cd _1 or one of its aliases. This frequency 
fulfilled the amplitude requirement in B, but not in V and 
only barely in U. It must therefore be regarded as a candi- 
date frequency only without further observational confirma- 
tion. Moreover, we could not distinguish the frequency from its 
aliases from multiperiodic fits to the lightcurves. Unlike for the 
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Fig. 1. The periodograms for the Geneva U (left) and B (right) data of HD 203664 after subsequent stages of prewhitening with 
the frequencies provided in the text. The amplitudes are expressed in magnitudes. The dashed horizontal line indicates the 4cr 
level determined in the way explained in the text. 
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Fig. 2. Phase diagrams of the U (left), B (central), V (right) lightcurves of HD 203664 for /i (upper) and for f 2 after prewhitening 
with /i (lower). The dots are the data and the full line is a sinusoidal fit with fixed indicated frequency. Note the different scale 
of the y-axes of the upper versus lower plots. 
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second frequency, numerous aliases of /a gave equally good 
fits. 

The best overall fits to the U, B, V data was achieved by 
taking the values for fx , . . . , / 3 as listed above, but we cannot 
exclude to have missed a true frequency and have taken an alias 
for any of f 2 , f^. We will not use the latter frequency to make 
seismic inferences further on and in the case of the second fre- 
quency we consider each time both f 2 and f 2 . 

Any other frequencies found after prewhitening do no 
longer reach four times the noise level (see bottom panels of 
Fig-Q so we stop the frequency analysis at this point. 

A firm conclusion is that HD 203664 is a multiperiodic 
(3 Cephei star with one dominant mode having an amplitude 
about ten times larger than the ones of its other modes. The 
three frequencies fx, f 2 , f% reach an amplitude of respectively 
38.8, 4.5, 4.0 times the noise level in U, of 38.4, 4.7, 4.2 times 
the noise level in B and of 31.8, 4.5, 3.2 times the noise level 
in V. 

4. Amplitude ratios and mode identification 

Table^lists the results for the amplitudes and phases of least- 
squares fits with fx , /g, /a fixed at their values mentioned 
above. The overall variance reduction is also listed. Using f 2 
rather than f 2 , or any alias of fa rather than these frequencies, 
led to amplitude and phase values within the error bars listed 
in Tabled As a comparison, we mention that a monoperiodic 
fit with fx has a variance reduction which is typically 5 to 8% 
lower depending on the filter. The amplitudes were each time 
largest in the U filter, as is expected for the oscillation modes 
in the (3 Cephei stars. 

We observed that the phases are equal to within their accu- 
racy for the three modes in the seven filters. Hence, we made 
use only of amplitude ratios as a mode identification diagnostic. 
This is according to the common procedure for the Cephei 
stars. 

In order to attempt mode identification, we proceeded as 
follows. We computed stellar models using the Code Liegeois 
d'Evolution Stellaire (CLES) kindly provided by R. Scuflaire. 
For the details on the input physics, used opacity tables and 
equation of state, we refer to Ausseloos et al. (2004). We re- 
stricted to models with X = 0.70, Z — 0.02, in agreement 
with the abundances derived by Little et al. (1994), and without 
core convective overshooting. For each value of the mass, we 
computed evolutionary tracks from the ZAMS until the TAMS 
and we selected those that are within the observed range of T e g 
and logg derived in Sect. 2. In doing so we considered a range 
in mass from 12 to 16 M® in steps of O.5M . Subsequently, we 
computed eigenfrequencies and eigenfunctions for each of the 
models using the non-adiabatic oscillation code MAD (Dupret 
2001) kindly made available by M.-A. Dupret. For each of 
the models, we selected the eigenfrequency which was clos- 
est to the measured fi, i — 1, . . . , 3 and considered its theo- 
retical non-adiabatic amplitude ratios following the formalism 
by Dupret et al. (2003). Finally, we compared all these theo- 
retically predicted amplitude ratios with the observed ones for 
each of the detected frequencies. In this way, we performed a 
mode identification which is independent of one particular stel- 



lar model, but that considers a large range of theoretical models 
covering safely the current mass estimate of HD 203664. 

The result of this procedure for the dominant mode can 
be found in Fig.[5]for £ = 0, . . . , 4. Thanks to the small er- 
ror bars of the observed ratios, we readily identified this mode 
as a dipole mode. Indeed, the I — 2 solution was the only one 
compatible with the high-quality data. The strongest mode of 
HD 203664 is therefore clearly non-radial. A similar situation 
occurs in several other (3 Cephei stars (see, e.g., Heynderickx 
et al. 1994). 

The other two modes had too high errors on their ampli- 
tude ratios to be able to discriminate among the ^-values. We 
illustrate this for the mode with frequency f 2 in FigTJ] While 
the agreement between the theoretical predictions and the ob- 
servations is best for £ = 2, we cannot firmly exclude the other 
^-values of the non-radial modes. We do find that the mode 
with frequency f 2 is unlikely to be radial. The same result is 
obtained if we take f 2 rather than f 2 . 

The ratios for the modes with f$ are even more uncertain 
and do not provide any constraint at all so the corresponding 
mode identification plot is omitted here. 

In an attempt to confront these seismic observational con- 
straints with standard models, we used the ones mentioned 
above, i.e. with masses ranging from 12 to 16 M , X = 
0.70, Z = 0.02 and without core overshooting (a ov = 0.0). 
We computed all their oscillation frequencies for zonal modes 
of £ = 0, . . . , 3 without taking into account the effects of ro- 
tation. We point out that this is a simplification of the true situ- 
ation, because second-order rotational effects imply frequency 
shifts, even for zonal modes. Moreover, mode coupling occurs 
due to rotation for modes whose degree differs by 2 if their fre- 
quencies are close together (e.g. Soufi et al. 1998, Daszynska- 
Daszkiewicz et al. 2003). However, we postpone a sophisti- 
cated interpretation of HD 203664's observed frequencies for 
the time being and restrict to a confrontation with standard 
models here. 

Typically three or four radial orders of an I = 2 zonal mode 
fitting the frequency fx resulted from the standard models from 
the ZAMS to the TAMS, for each of the considered masses be- 
tween 12 and 16 Mq. Forcing a simultaneous fit to an I — 2 
zonal mode for f 2 was achieved for several of the models that 
already fit fx- An examples of the mode spectrum of evolution- 
ary models with M = 13 M© is provided in Fig.|5] in which all 
the £ = 0, . . . , 3 mode frequencies are compared with fx and 
fi. The vertical dashed-dot line indicates the lower limit of the 
observed effective temperature interval, i.e. HD 203664 must 
occur to the left or close to this line, pointing out its young 
age already found by Aerts (2000). The model indicated by 
the full vertical line provides a good fit to fx and f 2 for zonal 
£ — 2 modes. This model has parameters X ~ 0.70, Z — 0.02, 
M = 13 M , log T off =4.429, log.g = 3.98, logL/L =4.24, 
R = 6.12 Rq, a central hydrogen abundance X c = 0.438, 
an age of 7.5 x 10 6 yr and no core overshooting. These T e g 
and log g values are compatible with the estimates from obser- 
vations provided in Sect. 2. Moreover, non-adiabatic computa- 
tions with MAD indicate that the two modes are excited for this 
model. 
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Table 1. Results of harmonic fits to the Geneva lightcurves of HD 203664. A stands for the amplitude, expressed in millimag 
and <j> for the phase, expressed in 2tt radians. The adopted reference epoch for <\> = 0.0 corresponds to HJD 2450000.0. 
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Fig. 3. Observed amplitude ratios Ax/Ajj (filled circles) and their uncertainties for the dominant mode with frequency /j of 
HD 203664, where Ax stands for any of the amplitudes in the seven Geneva filters U, Bi, B, B2, Vi, V, G. The lines represent 
theoretical predictions in the non-adiabatic treatment of the oscillations for stellar models within the mass range [12, 16] Mq 
during the main-sequence phase and fulfilling the error box in T c g and log g derived in Sect. 2. 
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Fig. 4. Same as Fig. [3] but for the mode with frequency fi. 
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Fitting /1 and f 2 as two zonal £ = 2 modes is also possible 
if their radial order differs by 2, but the models for which such 
a fit is obtained are somewhat more evolved which is less likely 
due to the observational constraint on the effective temperature. 
Moreover, the frequency fit is not as good as for Any of the 
models with M 7^ 13 M© while keeping X, Z, a ov fixed do 



not provide an equally good fit to both fx and f 2 (or f 2 ) as the 
model indicated by the vertical line in Fig.|5] 

We thus find that a good fit for fi and f 2 as I = 2 modes 
can only be achieved close to an avoided crossing if we restrict 
to zonal modes. Modes near avoided crossings have very good 
potential to probe the interior of a star. However, we have no 
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Fig. 5. Theoretical frequency spectra for stellar models with M — 13 M Q , X = 0.70, Z = 0.02 without core overshooting from 
the ZAMS until the TAMS. The star's position is to the left of the vertical dashed-dot line according to the T c s -range derived in 
Sect. 2. The symbol convention for the modes is as follows: squares for £ = 0, triangles for I = 1, plus sign for I = 2, circles for 
I = 3. Higher values are not shown for clarity. The observed frequency values f\ (£ = 2) and / 2 (likely I = 2) are indicated 
by the horizontal lines. The vertical line indicates the model for which an exact match to fi and / 2 is obtained for two zonal 
I — 2 modes. 



conclusive observational evidence to restrict to m — and so 
other evolutionary models may also lead to an equally good fit 
of the two considered frequencies for non-zonal modes. 

At present we have too restrictive observational constraints 
to perform a precise fitting procedure (e.g. as outlined in 
Ausseloos et al. 2004), i.e. to scan in detail the parameter space 
of any possible stellar model as a function of (X, Z, a ov , 
M). Indeed, besides the fact that we cannot exclude a differ- 
ent t— value for the mode with frequency / 2 (or / 2 ), we have 
no information on the azimuthal number of the modes nor do 
we know the surface rotation frequency of the star. Assuming 
that f\ and / 2 represent to a good approximation the central 
peaks of frequency multiplets is too restrictive because the ro- 
tational frequency must be of order 0.5 cd -1 or larger given 
that vsini=200kms^ 1 . 

We can only conclude that the currently observed oscilla- 
tion spectrum of HD 203664 can easily be explained by stan- 
dard stellar models and excitation computations and that we 
have insufficient information to further constrain the stellar pa- 
rameters from asteroseismology at present. 



5. Discussion 

We have shown the large-amplitude (3 Cephei star HD 203664 
to be a multiperiodic non-radial oscillator from single-site 
long-term multicolour Geneva photometry. The dominant 
mode was unambiguously identified as a dipole t = 2 mode 
from its amplitude ratios. The frequencies of the star are com- 
patible with standard stellar models of massive stars during the 
main sequence phase. 

There are at present nine other (3 Cephei stars known 
to have a peak-to-peak V amplitude larger than the one of 
HD 203664 (Stankov & Handler 2005). Only one of these ten 
stars is monoperiodic (BW Vul) since multiperiodicity was re- 
cently discovered in HD 180642 (Aerts et al., in preparation). 
Three of them have a confirmed dominant radial mode (v Eri - 
Handler et al. 2004, BW Vul - Aerts et al. 1995, HD 180642 - 
unpublished), three have an I = 1 mode (IL Vel - Handler et al. 
2003, 12 Lac - Handler et al. 2005, KPPer - Saesen et al., in 
preparation) and two have an i = 2 mode (KZ Mus - Handler 
et al. 2003, HD 203664 - this paper). We do not find any obvi- 
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ous relation between these results and the projected rotation ve- 
locity of the stars (Stankov & Handler 2005). HD 203664 is by 
far the most rapid rotator among them, with vsini=200kms~ 1 
which is more than twice as high as for any of the other large- 
amplitude (3 Cephei stars. From the mass and \ogg estimates 
given in Sect. 2, we find a critical velocity of 480 km s -1 and 
hence vsini / v cr ;t = 42%. 

The stellar parameters and internal structure of HD 203664 
can only be constrained further by means of elimination of the 
alias problems presented here for the low-amplitude modes (i.e. 
from a multisite campaign as those performed by Handler et al. 
2003, 2004, 2005), by subsequent unambiguous identification 
of their degree from amplitude ratios and by further identifica- 
tion of their azimuthal number through high-resolution spec- 
troscopy. Given its visual magnitude of 8.6, its oscillation pe- 
riods near 4 hours and its low-amplitude modes, together with 
the data requirements for such type of spectroscopic analysis 
(e.g. Aerts & Eyer 2000), one needs, besides a photometric 
multisite campaign, at least one week of 8-m class telescope 
time to achieve mode identification and subsequent interpreta- 
tion of the oscillations of this massive hot rapid rotator among 
the j3 Cephei stars. 
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